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ABSTRACT. We consider the Hartree equation with a smooth kernel and an external po- 
tential, in the semiclassical regime. We analyze the propagation of two initial wave packets, 
and show different possible effects of the interaction, according to the size of the nonlin- 
earity in terms of the semiclassical parameter. We show three different sorts of nonlinear 
phenomena. In each case, the structure of the wave as a sum of two coherent states is pre- 
served. However, the envelope and the center (in phase space) of these two wave packets 
are affected by nonlinear interferences, which are described precisely. 
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1. Introduction 
Consider the following Hartree equation in the semiclassical regime e — » 0: 

£ 2 

(1.1) i£S# + -Af = V(t,x)ijf + e a (K* \i/j E \ 2 )ip e , teR + = [0,oo), x <= TL d , 

where a > 0, K : H d -> R, V : R+ x R d -> R, d ^ 1. Equation CD appears for 
instance as a model to study superfluids, with application to Bose-Einstein condensation: 
in 1 5 6 1, the kernel K is given by the formula 

K(x) = («i +a 2 \x\ 2 + a 3 |x| 4 ) e^^ 1 ^ + a 4 e~ s2|:r|2 , a u a 2) a 3 , a 4 , A, B e R. 
Assume 

(1.2) V £ (0,x) =e- d,A a (^j=^J e l{x -" o) - po/E , a £ S{R d ), qo,po G R d . 

Such initial data are called semiclassical wave packets, or coherent states. They correspond 
to a wave function which is equally localized in space and in frequency (at scale \/e), so 
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the uncertainty principle is optimized in terms of e: the three quantities 

II^(o)||l>(r<), 



, and 

L 2 (R d ) 



— 7=-^ (0) 



L 2 (R d ) 



have the same order of magnitude, 0(1), as e — > 0. In the linear case K = 0, another 
reason why such specific initial data are particularly interesting is that the superposition 
principle is available: if we can describe tp s in the case (11.2b . then the evolution of a sum 
of initial wave packets of the form ( 11.2b is simply the sum of the evolutions of each initial 
wave packet. In this paper, we address this question in the nonlinear setting. We describe 
several nonlinear interference phenomena in the case where K ^ is smooth, and ip E (0, x) 
is the sum of two such wave packets. 

The value of the parameter a in ( II . lb measures the strength of the nonlinear interaction 
in the limit e — > 0. In [ 12 1, where the Hartree nonlinearity is replaced by a local nonlinear- 
ity, it is established that if nonlinear effects are critical in terms of semiclassical dynamics 
(that is, the value of a is critical, see §1.2\ for this notion), then despite the fact that the 
problem is nonlinear, the superposition principle remains valid, in the limit e — > 0. In |9), 
the case of a homogeneous Hartree nonlinearity K(x) = A|a;|~ 7 is considered: conclu- 
sions similar to those in STU are proven. In these two frameworks, the description of the 
wave packet dynamics in a "supercritical" case (nonlinear effects are stronger than in the 
critical case) is an open question, even on a formal level. On the other hand, in the case 
of a smooth Hartree kernel, the propagation of a single wave packet has been described in 
supercritical regimes (|]3]|9]). In this paper, we prove that in the critical regime, nonlinear 
interferences affect the propagation of two initial wave packets at leading order, in contrast 
with the case of a homogeneous kernel. We also describe the nonlinear interactions in su- 
percritical regimes, where even stronger interferences are present. In all cases, we prove 
a convergence result on all finite time intervals (t € [0, T] with T independent of e), as 

£ -> 0. 

Assumption 1.1. The external potential V is C 3 , real-valued, and at most quadratic in 
space: 

V E C* 3 (R+ x R d ; R), and d?V e L°° (R+ x R d ) , \f3\ = 2, 3. 
The kernel K is C 3 , real-valued, bounded as well as its first three derivatives: 

K e C* 3 n VK 3 <°°(R d ;R). 
Consider the Hamiltonian flow: 

(1.3) q(t)=P(t), p(t) = -VV(t,q(t)); q(0) = q , p(0) = Po . 
The regularity of V implies that (11.3b has a unique, global solution 

t-> (q(t),p(t)) e C 3 (R+;R 2d ). 

Since we shall consider only bounded time intervals in this paper, the growth in time of the 
classical trajectories is not discussed. 

1.1. The linear case K = 0. Introduce the function 

cpUt,x) = e-^v}™ (t, ^§^) e ^(*)+p(*)-(*-<K*)))/-, 

where (q,p) is given by ( 11.3b , the classical action is given by 

(1.4) S(t) = J Qb(s)| 2 - V(s, q(a))j ds, 
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and the envelope u = u (t, y) solves 

(1.5) id t u Ym + l -kv}™ = ±(y,V 2 V(t,q(t))y)u Vla ; u lia (0, y) = a{y), 

where the notation V 2 stands for the Hessian matrix, and since the space variable for u hn 
is y, A stands for A y . The following lemma is standard, see e.g. ll4l [l3l[T4l[T5ll24ll25ll26l 
and references therein. 

Lemma 1.2. Let a e 5(R d ), and tp e solve ( 11.11 ) with K — 0, and ( II. 21 ). There exist 
positive constants C and C\ independent of e such that 

W{t) - ^ in (t)\\ LHRd) < C^e c ^, Vt > 0. 

In particular, there exists c > independent of e such that 

sup \\^ 6 {t) - ^f in (f)|| L 2 (R£i) — >0. 

1.2. Nonlinear case: notion of criticality. In the nonlinear case K ^ 0, the following 
distinction was established in J9J: 

• If a > 1, nonlinear effects are negligible at leading order: with the same function 
<y9f in as in the previous section, there exists C > such that 



v 



If a = 1, nonlinear effects become relevant at leading order (unless K(0) = 0): 
there exists C > such that 

L 2 (R d ) 

From this point of view, the case a = 1 is critical: the supercritical behavior is described 
in two cases, a = 1/2 (|9|) and a = ([3 9 |). The approximate solution derived in these 
two cases may be viewed as a particular case of the approximate solution presented below, 
when one of the two initial wave packets is zero, so we choose not to be more explicit 
about these two cases here. Other cases could be described as well: the case a e (0, 1/2) 
is similar to the case a — 0, and the case a <G (1/2, 1) is similar to the case a = 1/2, up 
to several modifications in the notations essentially. 

In the case a > 1, nonlinear effects are negligible at leading order, so the superposition 
principle remains: the nonlinear evolution of two (or more) initial wave packets is well 
approximated by the sum of the linear evolutions of each wave packet. We will see that 
when a ^ 1, nonlinear interferences affect the behavior of ip e at leading order. 

Throughout this paper, for k E N, we will denote by 

S fe = | / e L 2 (H d ) ; ||/|| E > := £ ||^/|| i2(Rd) < oo 

and E 1 = S. As established in |9), if ip E (0, •) € L 2 (R d ), then under Assumption [TT] 
( II. U has a unique solution ip E £ C(R+; L 2 (R d )), regardless of the value of cfl and 

||V> e (*)l|L»(R-) = 11^(0)11^(114), vt^o. 



'To be complete, the regularity assumption on V in (9) is stronger, but Assumption |l.l| is enough. 
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1.3. Critical case: a = 1. We now consider ( 11.1b in the case of two initial wave packets: 
(11.2b is replaced by 



(1.6) 



V> £ (0,x) 



3=1,2 



with (gio,pio) 7^ (920, P2o)- Let (oj . Pj) be the solution to (11.3b with initial data (q j0 ,Pjo) 
and the associated classical action given by (11.41 ). Define the approximate solution as 



(1.7) 



3=1,2 



<&■(*) 



^(Sj(t)+Pi(t)-(^-9j(t)))/£ 



where the envelopes Uj are given by the formulas: 

J 2/) = -^i in (t, 2/ ) e -^^C°>ll-illi 2 -' i ll-2ll^ 2 J r f J p 4 r( (?1 ( s )-^ 2 c S ))^ 7 
\u 2 (t,y) = 4 in (t,t / ) e - itiC ( )ll^lli2- i ll a ill^/ ' x feW-9iW)^ ! 
with obvious notations adapted from ( 11.5b . 



(1-8) 



Theorem 1.3. Lef d ^ 1, V, K satisfying Assumption \l.l\ Let a\,a% £ S 3 , and ip £ be the 
solution to (II. lb vv/f/; a = 1 and initial data (11.6b . Then for any T > independent of e, 
there exists C > independent of e such that 



te[o,T] 



where ipl pp fry ( 11.7b -( [TT8"l l. 



The nonlinear effects are described at leading order by the exponentials in d 1 . 8b - Even in 
the case of a single initial wave packet (say a 2 = 0), the nonlinearity affects the envelope 
by a phase self-modulation. The second terms in the exponentials describe the effect of 
nonlinear interferences, which are not a simple superposition in general. 

As pointed out above, it may be surprising to notice that even in the critical case a = 1, 
nonlinear interferences are present at leading order. This is in sharp contrast with the case 
of an homogeneous kernel, K(x) = A|x|~ 7 , < 7 < min(2,d). It was shown in (9) 
that in this case, the critical value for a is a c = 1 + 7/2, and that when a = a c , the 
superposition principle remains, even though the nonlinearity affects the propagation of a 
single wave packet at leading order (the envelope equation is nonlinear). 

1.4. Case a — 1/2. The approximate solution is now constructed as follows. The pairs 
(qj,Pj), j = 1, 2, are still given by the usual classical flow ( 11.3b . On the other hand, we 
modify the actions, and make them e-dependent: 

rt 



(1-9) 



sm 



l\ Pl {s)^-V{sMs)))^ 



tVi^(0)||ai||| 2(Rd 
1 



V^||02||i2(Rd) 

V(s,q 2 (s))] ds 



t^K(0)\\a 2 \\ 2 LH 



R d) - \/e||ai||i2(Rd) / 
Jo 



K( qi ( S )-q 2 ( S ))ds, 



K (^(s) - qi{s)) ds. 
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Consider the system of Schrodinger equations 



(1.10) 



idtui + ^Au x = - (y, V 2 y (t, qi (t)) y) u x + |K||| 2 y • Vtf(0)i2i 



1 



+ \\a 2 \\ 2 L2 yVK{q l {t)-q 2 {t))u ll 
1 



id t u 2 + ^Au 2 = - (y, V 2 V (t, q 2 {t)) y) u 2 + \\a 2 \\ 2 L2 y ■ VK(0)u 2 

+ \\a 1 \\l 2 yVK(q 2 (t)-q 1 (t))u 2 , 

with initial data ai and a 2 , respectively. These are two linear equations with time depen- 
dent potentials, which are polynomial in y, of degree (at most) two. The following result 
is classical, see e.g. ll36l[T7l[T8l : 

Lemma 1.4. For j = 1,2, let a 3 - e L 2 (R d ), and (q 3 ,p 3 ) € G 3 (R + ;R 2d ) given by 
< l 1 - 3b . There exists a unique solution (ui,u 2 ) £ G(R+; L 2 (R d )) 2 to ( II. 10b such that 
(ui,W2)|t=o = (oii a 2)- In addition, the following conservations hold: 

\\uj(t) \\l'(b.*) = ||a J || i 2( Rti) . Vi ^ 0, j = 1,2. 



To define the envelopes in ( 11.7b . set 
(1.11) 



u 2 (t,y 2 ) = u 2 (t,y 2 )eyip i 



t»i(t, yi) = fii(t, yi) exp ( t / [Vif(0) • G x (s) + Vif (gi(s) - g 2 (s)) • G 2 (s)J dsj , 

(Vif(0) • G 2 (s) + Vif (q 2 (s) - qi {s)) • Gi(s)) ds^j , 

where = / z\iij(t, z)\ 2 dz. Since Gj is a nonlinear function of Uj, the system 

formed by (1*1,142) is nonlinear, with a nonlinear coupling: nonlinear interferences are 
present both in rapid oscillations — the modified actions generate y^-oscillations in time 
— and in the envelopes. The presence of the functions Gj in the above formulas reveals 
non-local (in space) nonlinear phenomena concerning the envelopes in tp^ pp . Since the 
problem is now supercritical, it should not be surprising that stronger regularity properties 
are assumed in the following result (see Remark |4~2V 

Theorem 1.5. Let d 1. Assume that V and K are real-valued and satisfy: 

V E G 5 (R + x R d ;R), and d^V £ L°° (R + x R d ) , 2 < |j8| < 5. 
K e W^ 5 '°°(R d ;R). 

Let ax,a 2 G S 5 , and ^ pp be given by (IQ-rfL^- CTOl l- CTB . Then for any T > 
independent of e, there exists G > independent of e such that 

sup \\r{t)-^l w {t)\\ <GVi. 
te[o,T] ^ ; 

1.5. Case a = 0. In this last case, nonlinear interferences affect even the geometric prop- 
erties of the wave packets, in contrast with the cases a = 1 and a = 1/2. The trajectories 
are required to evolve according to the system 

' qi(t) =pi(t), 

Pl (t) = -VV(t,q 1 (t))-\\a 1 \\ 2 L2 VK(0) 
hit) =p 2 (t), 
{ p 2 (t) = -W (t, q 2 (t)) - ini^vi^o) 



(1.12) 



\a 2 \\ 2 L2 VK ( qi (t) - q 2 (t)) . 
\ ai \\ 2 L2 VK(q 2 (t)- qi (t)). 
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Unless VK is a constant (which would implies that K is constant, a trivial case), one can- 
not decouple the unknowns (qi,pi) and {q 2l p 2 ): the coupling cannot by undone, and the 
"good unknown" is (qi,px, q 2 ,p 2 ) € R 4d . In view of Assumption ll.il Cauchy-Lipschitz 
Theorem implies: 



Lemma 1.6. For j = 1,2, let (qjo,pjo) S R 2d . If V and K satisfy Assumption U.l] then 
(11.12b has a unique solution (<jfi , pi , q 2 ,p 2 ) € C 3 (R+;R 4d ). 

Remark 1.7 (Hamiltonian structure). If the external potential does not depend on time, 
dtV = 0, and the kernel K is even, K(—x) = for all x £ R d , then the Hartree 

equation (II . lb has a Hamiltonian structure. In the case a = 0, the following energy is 
independent of t, 

Note that since -ftf is even, Vif (0) = 0, and VK{q 2 - qi) = -Vif(gi - <? 2 )- In that 
case, the system of modified trajectories (11.12) is also Hamiltonian, as can be seen from 
the approach presented in [21 1. Given the state variable z — {q\,p\, q 2 ,p 2 ) T , l e t 



H(t,z) = ai -\ Pl \ 2 +V(qi) ) + a 2 [ -\p 2 \ 2 + V(q 2 ) + a x a 2 K {q x - q 2 ) , 



where a. 



*3 II L 2 



The system (11.12b has the Hamiltonian structure 



dz 



JD z H{t, z) with J 



( 





l/ai 





° \ 




-l/ai 






















l/«2 


\ 








-l/a-2 


/ 



One checks indeed that the function H is conserved along solutions of ( 11.12b . 

Before defining the modified actions, we have to construct also the envelopes. Consider 
the coupled, nonlinear system 



id tUl + ^A Ul = l - (y, M x {t)y) u x - (V 2 X(0)C7 1 (t), y) u x 

- (\7 2 K( qi - q 2 )G 2 (t),y) Ul + \ ( f (z,\7 2 K(0)z) \ Ul {t,z)\ 2 dz ) u x 



(1.13) { 



+ \(J ( z > v2 ^(9i - M*, z)\ 2 dz^j ui, 
id t u 2 + ~Au 2 = ~ (y, M 2 {t)y) u 2 - (V 2 K(0)G 2 (t), y) u 2 



(V 2 K{q 2 - q 1 )G 1 {t),y)u 2 + -U (z,V 2 K(0)z) \u 2 (t,z)\ 2 dz ) u 2 
X - (J (z,W 2 K(q 2 - qi )z) \ Ul (t, z)\ 2 dz^j u 2 , 
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where the functions qj are assessed at time t, and we have denoted 



(1.14) 
(1.15) 



G 3 (t) 
Mi(t) 



R'' 



z\uj(t,z)\ 2 dz, i = l,2, 



a illi 2 (R^)V 2 ^(0) + ||a 2 || 2 L2(Rti) V 2 if ( gi (i)- 92 (i)) 
+ V^(t,«i(t)), 

(1.16) M 2 (t) = ||a 2 |li 2(Rti) V 2 X(0) + ||a 1 || 2 L2(Rd) V 2 ^( g2 (i)- 9l (i)) 

The system defining the envelopes is more nonlinear than the cases a = 1 and a = 1/2, 
and, as in the case a = 1/2, involves nonlinear terms which are non-local in space. In 
Section|5] we prove the following result: 

Proposition 1.8. Let {q\ , p\ , q 2 , p 2 ) be given by Lemma \LE\ and a\ , a 2 € S fc w/f/i fc 1. 
Then ( I1.13l l /ias a unique solution (ui,it 2 ) G C(R+;X fc ) vWf/i initial data (ai,a 2 ). /n 
addition, the following conservations hold: 

||%(*)IU 2 (Rrf) = ll a jlU 2 (R d )< Vt ^ 0, j = 1,2. 
We can then define the modified, e-dependent actions: 



[ (^Ms)l 2 - V (s, qi (s)) - #(0)110x111* - K( qi (s) - q 2 (s)) \\a 2 \\ 2 L2 
+ y/2VK(0) ■ Gx{s) + V~eVK {q x {s) - q 2 {s)) ■ G 2 (aj)da, 
(\\P2(s)\ 2 - V (s,q 2 (s)) - ^(0)||a 2 || 2 2 - K (q 2 (s) ~ qx(s)) \\ai\\ 2 L2 
+ ^VK (0) • G 2 (s) + ^feVK (q 2 {s) - q x {s)) • d(«)W 



Sf(t) 



Theorem 1.9. Let d ^ 1. Assume that V and K are real-valued and satisfy: 

V G C 7 (R+ x R d ; R), and d?V G L°° (R+ x R d ) , 2 < < 7. 
# G VK 7 '°°(R d ;R). 

Lef ax, a 2 G £ 7 . There exist 9\, 6 2 G C 2 (R+; R), w;f/i 0j(0) = 0j(0) = 0, such that the 
following holds. For any T > independent of e, there exists C > independent of s such 
that 



sup 

te[o,T] 



3=1,2 



L 2 (R d ) 



where we have denoted 



^.(t,x) = e~ d l\ (t, X ~J W ^) e <s?W+P,(tM*-*M))/e 



In general, the phase shifts 6*j are not identically zero. In Section|S] we give the expres- 
sion of these functions, which is probably a bit too involved to present at this stage (see 
Equation ( 18.4b ). and check that in general, (#i(0), # 2 (0)) ^ (0, 0). Such modifications do 
not appear in the case of a single wave packet studied in J3]|9]. Technically, the reason is 
two-fold. First, in EH], it is assumed that WK(0) = 0, so the first line in (18.4b vanishes. 
Then, the second line in (18.4b accounts for the presence of two wave packets, and measures 
some coupling through a phase modulation: it vanishes in the case of a single wave packet. 
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1.6. Comments. 

The results. In the three cases studied here, the interferences are nonlinear. They always 
affect the envelopes. In the case a = 1/2, they affect moreover the action, and in the case 
a = 0, the system ( 11.12b reveals a nonlinear coupling of the trajectories, so all the terms 
involved in V'app are influenced by the nonlinearity. Our results are original even in the 
case V = 0. 

Nonlinear interferences always carry a non-local in time aspect. Even if K is decaying 
at infinity, the interactions ignore the mutual distance of the two wave packets: no matter 
how large qi (t ) — q-i (t) is, nonlinear interferences affect the solution at order O ( 1 ) on finite 
time intervals, as e — > 0. 

Our results yield a unified picture concerning Wigner measures (see e.g. |8 20 32]): 

Corollary 1.10. In all the cases a. = 1, a = 1/2 and a = 0, and under the Assumptions 
of Theorem U 3\ Theorem U .5\ and Theorem U .91 respectively, the Wigner measure ofip e is 
given by 

w(t,x,(,)= hj Whin*) 6 ( x -&(*)) ®<K£ ~Pj(*))> 

3=1,2 

with (qj , pj ) given by the standard Hamiltonian flow ( 11. 3t in the cases a = 1 and a = 1/2, 
and (qi,Pi,q2,P2) given by dl.l2t in the case a = 0. 

To be complete, the proof of this corollary relies also on the results established in Sec- 
tion [3] In the two cases a = 1 and a = 1/2, the Wigner measure of iptpp i s not affected 
by the nonlinearity, even though we have seen that the Hartree nonlinearity does affect the 
leading order behavior of the wave function, and that nonlinear exchanges are present at 
leading order. In the case a = 0, nonlinear effects alter the Wigner measure, even when 
VK(0) = 0, a case which is often encountered in Physics (typically, K(—x) = K(x), so 
the Hartree nonlinearity has an Hamiltonian structure). In other words, the Wigner mea- 
sure of tp € is always affected by nonlinear interferences. This is in contrast with the case 
of a single initial wave packet, where the trajectory (q,p) is modified as if an electric field 
||a||| 2 VK(0) ■ x was added to the initial Hamiltonian -^A + V: if Vif (0) = 0, then the 
Wigner measure ignores nonlinear effect even in the case a = (see |[3][9]). 

Note that if Vif(0) ^ (a case which is not necessarily physically relevant, from the 
above remark), Theorem U. 9 l is new even in the case of a single wave packet. 

In this paper, we treat the case of two initial wave packets: our approach can easily be 
generalized to any (finite) number of initial wave packets, the main difference being that 
formulas get more and more involved as the number of initial wave packets increases (but 
the main analytical aspects are essentially the same). 

We have examined the leading order behavior of the exact solution, up to an error of 
order 0(y/e): like in Q, ip £ could be approximated by a series involving powers of ^Je, 
up to arbitrary order. This statement is made more precise in fJH] (see in particular Equa- 
tion ( 18.2b ): to prove Theorem |1.9l the asymptotic expansion of the main unknown functions 
has to be pushed one step further than in the cases a = 1 and a = 1/2. 

Comparison with related works. We briefly give more details concerning the propagation 
of two wave packets described in 1(9] 1121 . Since both cases are rather similar, we describe 
the case of a Hartree nonlinearity only ([9 1). The main difference in the setting is that dl.lt 
is replaced with 

ied t ip £ + —Aip e = V(t, x)il> e + Xe a (|a;|- 7 * \^ E \ 2 ) ^ e , t> 0, x e R d , 
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with A G R and < 7 < min(2, d). The critical value for a is a c = 1 + 7/2 > 1. When 
a = a c , the propagation of one initial wave packet is well approximated by 

<f(t,x) = e-^u (t, e *(SW+P(t)-(«-9(*)))/ e) 

where (q,p) is given by dl.3t , S is the classical action dl.4K and the envelope u solves the 
nonlinear equation 

id t u + ^Au = i (y, V 2 V (t, q(t)) y)u + X flyp * \u\ 2 ) u. 

It is proved that two such wave packets evolve independently from each other, up to an error 

which is O 2(1 +t> ^ . A way to understand this result compared to the ones presented here 

is that since a c > 1, no interference can occur at leading order. 

There are several results which may seem similar to ours, in the case of one initial wave 
packet: see e.g. Q [16] [19] [27] [28] [29] [30]. In those papers, the initial amplitude a is 
very specific, since it is a ground state. The propagation and stability of multi-solitons 
for the nonlinear Schrodinger equation (without external potential) have been studied in 
Il33l[34l[35ll37l (see also ll40l ). In the framework of these papers, the waves do not interfere. 

In 0]], a problem which shares several features with ours is studied: there is an external 
potential, the regime is semiclassical (see [27 1), and nonlinear. The envelopes of the initial 
data are two soli tons. The structure of the soli ton manifold implies some rigidity on the 
evolution of the initial data. Eventually, the two waves do not interact at leading order. 

On the other hand, in (3T), the case of two soli tons for the Hartree equation has been 
studied. In this non-semiclassical setting, and in the absence of an external potential, the 
authors construct a solution which behaves, for large time, like the sum of two solitary 
waves, whose respective centers in phase space evolve according to the two-body problem. 
This feature can be compared to Theorem ll ,9l (with V = 0), where the centers of the wave 
packets evolve according to the nonlinear system ( 11.12) . Nevertheless, the envelopes are 
given by the ground state, and do not evolve with time. The analytical approach is different: 
in (3Tj, a fine study of the Hartree operator linearized about the soli ton is performed, in 
particular to understand the spectral properties of this operator. On the other hand, we 
do not consider such an operator; a similar approach with general profiles 01, a% like we 
consider would probably be out of reach. 

In ||23l[38l[39l , a semiclassical regime is studied, in the presence of an external potential 
and a nonlinearity. The potential is a double well potential, and the associated Hamiltonian 
has two eigenfunctions. For initial data carried by these two eigenfunctions, it is shown 
that the nonlinear solution may remain concentrated on the eigenfunctions, with time- 
dependent coefficients which interact nonlinearly. 

In all the cases mentioned above, the nonlinear interference of the envelopes is negligi- 
ble, due to the fact that the envelopes decay exponentially. In our case, the decay may be 
much weaker (algebraic). However, even though we have seen that the envelopes always 
interact nonlinearly in the cases studied here, we will see that some "rectangle" terms are 
negligible in the limit e — > 0, thanks to a microlocal argument (see Section[5]l. 

We finally point out that nonlinear interactions of amplitudes have been analyzed in the 
context of weakly nonlinear geometric optics for Schrodinger or Hartree equations in var- 
ious contexts (not to mention the even wider literature concerning hyperbolic equations); 
see for instance ifTTl l22l . 
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Notations. Throughout this paper, R + stands for [0, oo). We also use the standard con- 
vention, for A 6 R™, n ^ 1, 

(A) = VTTW- 

For two positive numbers a £ and b £ , the notation a £ < b £ means that there exists C > 
independent of e such that for all e e (0, 1], a £ < C6 E . 

2. Formal derivation 

We resume the same approach as in the case of a single wave packet ((H), in the case 
of ( 11.61 ): from this point of view the computations below include the ones presented in J9j- 

2.1. The general strategy. We seek an approximate solution of the form 
(2.1) iftwfrx) = e" d/4 £ Uj (t, X ~ q ' {t) ) e *(^W+PiW-(*-*W))/=, 

for some profiles Uj independent of e, and some functions Sj (t) to be determined. These 
functions Sj correspond to the classical action (11.4) in the linear case. We will see that 
according to the value a in Jl.lt . the expression of Sj may vary, accounting for nonlinear 
effects due to the presence of the Hartree nonlinearity, and so it may be convenient to 
consider e-dependent functions Sj. Also, according to the value of a, the pairs (qj,Pj) 
will solve the standard Hamiltonian system (11. 3t , or a modified one. Denote 

<t>j{t,x) = Sj{t)+ P j{t)-{x-qj{t)). 

In the cases a = 0, 1/2 and a = 1 considered in this paper, we will see that we can write 

feftl^pp + yAVIpp - lAAapp ~ £ a (K * |^ pp | 2 ) Vfpp - 

(2.2) 



-d/4 



: ■ ■ 2, e*^/ £ (6 0j + ^i,+e6 2j +er|) (t, - ^ 
j=i, 2 V v£ 



for 6y independent of e. The approximate solution ipl pp is determined by the conditions 

boj = b hj = b 2 j = 0, j = 1,2. 

The remaining factor rj accounts for the error between the exact solution V> £ and the ap- 
proximate solution ip £ pp . Introduce two new space variables, which are naturally associ- 
ated to each of the two approximating wave packets: 

y^ x -jM, ; = i,2. 

At this stage, the introduction of these variables may seem very artificial, since only the x 
variable will eventually remain. It can be understood as a change of variable corresponding 
to the moving frame of each wave packet. Technically, it will be justified by the fact, 
already present in the linear case K = 0, that the remainders r| will satisfy pointwise 
estimates of the form 



t x-qj{t) 



<V~e(y 3 ) 3 A £ (t,y 3 ) 



Uj 



_ s-gj(t) ' 



The functions A £ - are well localized, in the sense that yj M- (yj) k A £ (t, y.j) in bounded in 
L 2 (R d ) at least for k = 3 (but possibly for larger fc's), while typically, a function of the 
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form 

cannot be controlled in £ 2 (R d ) uniformly in e and t £ [0, T] if rj ^ 0. 

To conclude this subsection, we expand each term on the left hand side of ( 12.21 i so it 
has the form of the right hand side. In the following subsections, we discuss the outcome 
according to the value a = 0,l/2ora = l. 

The linear terms are computed as follows: 



3=1,2 

a — 1 9 V / 



3=1,2 

Here, as well as below, one should remember that the functions are assessed as in ( 12. 11 1. 
Recalling that the relevant space variable for Uj is y 3 , we have: 

Mi = S 3 {t) + - ( Pj (t) ■ (x - q 3 (t))) = Sj(t) + V^Pjit) ■ Vj - Pj(t) ■ &(*)■ 
For the linear potential term, we write 

VF app = V(t,x)e- d '* Yl e* i{t ' x)fe Ui&Vi) 

3=1,2 

= e -«J/4 £ e ^-(M)/ £y ( M . (f)+y .^) 

3=1,2 

and we perform a Taylor expansion for V about x = q 2 (t): 
V (t, q 3 {t) + y jy /e) Uj(t, Vj) = V (t, q 3 {t)) u 3 \t, yj) + y/eyj ■ W (t, qj{t)) Uj (t, Vj ) 

+ \ (Vi , V 2 V (t , qj (t) ) Vj ) uj (t, yj ) + e 3/2 r e jV (t, Vj ) , 



with 



(2.3) \fjv&vj)\ <c< w ) 8 k-(*,yi)l, 

for some C independent of e, i and w,-, in view of Assumption ll.il In the case K = 0, we 
come up with the relations: 

b^J = -uj (sj(t) Pj {t) ■ qj(t) + + y ( t| 9i (t 

&y = -* (&(*) - Pi(*)) ' Vuj - % • (fe(t) + W (t, <&(*))) u.j. 



1 . 1 



6^ n = i9 tUj + — Awj - - (y jt V 2 V (t, qj(t)) y 3 ) u y 
For the nonlinear term, the computations are heavier: 

(K *m pp \ 2 )^l PP ^e- d ' A ]T e^ (t '* )/e ( / K(zm pp (t,x-z)\ 2 dz) Uj (t, yj ). 

3=1,2 VJ ' 

Eventually, each envelope Uj will solve a Schrodinger equation, the two equations being 
coupled. The precise expression of these equations depends on a, but at this stage, we can 
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notice that for j = 1,2, Uj solves an equation of the form 

(2.4) id t Uj + -Auj = - (yj , V 2 V (t, Qj (t)) y,) Uj + F jUj , 

where the function Fj, accounting for nonlinear effects due to the Hartree kernel, is real- 
valued. We infer an important property: the L 2 -norm of Uj is independent of time, 

(2.5) ll«i(*)lli»(R<«) = IKIU 2 (R.<*)' Vt^O, j = l,2. 
At this stage, this is only a formal remark. 

In the above sum, the variable x must be expressed in terms of yj : 



#*|V>app| 2 = / K(z)\^l pp (t,q 3 (t) + V^y 3 -z)\ 2 dz 



e- d ' 2 / K{z) 



J2 e iMt,*-*)/e Uk U y . + 
k=1.2 ^ 



Qj( t) - Qk(t) _ z_ 



dz. 



Before changing the integration variable, we develop the squared modulus: 



v e w-*y*u k ( t, Vj + gj(t) 7 gfc(t) - 4 

k=i,2 V Ve V£ 



ui t, y-j + 



qj(t)-qi(t) z 



u 2 t, yj + 



Qj(t) - q 2 (t) _ _z_ 



qj(t)-q 2 (t) z 



+ 2Ree^-W Ul (t, yj + - 2l) u 2 [t, Vj + 

where tfii — cf>2 stands for 4>i(t, x — z) — (f> 2 (t, x — z). To ease notations, we shall denote 
in the rest of this paper: 

Sq(t) = qi (t) - q 2 (t); Sp(t) = Pl (t) - p 2 {t). 

We can write 



3=1.2 



with 

V?\t, yi )=e- d l* lK{z)\ 



«i(^i-^ 



Sq(t) z 



2lW«- ^u 1 [t,y 1 -^)u 2 [t,y 1 + S M-^\ \d:. 



V 2 ^{t,y 2 )=e 



-d/2 



K(z) 



ui t, y 2 



+ 2Ree i( * 1 "* 2)/e ui [t,y 2 



Sq(t) z 

Sq(t) 



U2 t, y 2 



re ^ 



-7= )u 2 [t, y 2 7 =) \dz. 



Each nonlinear potential Vj is the sum of three terms. The third term in each of these two 
expressions, involving the product u\u 2 , will be referred to as rectangle term, as opposed 
to squared terms, involving squared moduli. The two rectangle terms are examined in 
Section [3] where we show that at least on finite time intervals, they are negligible in the 
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limit e — > 0, regardless of the value of a. Therefore, we now consider only the squared 
terms. Changing variables in the integrations and performing a Taylor expansion of the 
kernel K, we find successively (recall that Gj is defined by ( 11.14b ): 



-d/2 



K{z) 



ui [t,yi- 



dz — J K (y/e(yi - z)) \u\ (t, z)\ 2 dz 
K(0)\\aif L2 + V~e\\ai\\lm • VJC(O) - ^V^(O) ■ G x {t) 



+ -< yi ,V 2 K(0)j;i)|| a i||! 2 + - / (z,V 2 K(0)z)\ Ul (t,z)\ 2 dz 
-e(V 2 K{0)G 1 {t),y 1 )+e^ 2 f r e n (t, z - yi )\ Ul (t, z)\ 2 dz, 



-d/2 



K{z) 



u 2 [t,yi + 



Sq(t) 



dz 



K(5q(t) + V^(yi-z))\u 2 (t,z)\ 2 dz 
= K(Sq)\\a 2 \\ 2 L2 + Vi\\a 2 \\ 2 Lm ■ VK(6q) - J~eVK(5q) • G 2 (t) 
+ e -(y 1 ,V 2 K(Sq)y 1 )\\a 2 \\ 2 L2 + e - J (z,V 2 K{8q)z) \u 2 (t, z)\ 2 dz 

-e(V 2 K(Sq)G 2 (t), yi )+e 3/2 f r\ 2 {t, z - yi )\u 2 (t, z)\ 2 dz, 



-d/2 



K{z) 



Ml t, J/2 



Sq(t) 



dz 



- A-(-^)||oi||| a + ^IKHl^ • VK(-Sq) - ^IVK(-Sq) ■ Gi{t) 
+ i(y 2l V 2 K(-6q)y 2 )\\a 1 \\ 2 L2 + £ - J (z 1 V 2 K{~5q)z)\u 1 (t,z)\ 2 dz 

-e(V 2 K(-Sq)G 1 (t),y 2 )+e i/a / r s 21 (t, z - y 2 )\u 1 (t i z)\ 2 dz, 



-d/2 



K{z) 



U2 t, y 2 



dz 



= K(0)\\a 2 \\ 2 L2 + V~e\\a 2 \\ 2 L2 y 2 ■ VK(0) - y/eVK(p) ■ G 2 (t) 
+ i(y 2 ,V 2 K(0)y 2 )\\a 2 \\ 2 L2 + ^ J (z,V 2 K(0)z) \u 2 (t, z)\ 2 dz 

- e (V 2 K(0)G 2 (t), y 2 ) + e 3 ' 2 J r| 2 (t, z ~ y 2 )\u 2 (t, z)\ 2 dz, 
where the functions r| fc satisfy uniform estimates of the form 
(2.6) \r%{t,z)\^C{T){z)\ Vz e R d , t £ [0, T], 

with C(T) independent of e, j and k, but possibly depending on T. 
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2.2. The critical case: a = 1. When a > 1, we have bij = fe^ n for all £,j: nonlinear 
effects are not present at leading order. When a — 1, we still have bgj — bff for I = 0, 1: 
we still demand (qj,Pj) to solve ( 11.31 ) in order for the equations 6<y = 0, £ = 0, 1, to be 
satisfied, and Sj is defined as in ( 11.4b . On the other hand, the expression for b 2 j is altered: 

&21 = id t u! + -A«i - i (j/i, V 2 T^(t,qi(t)) yi) mi - A~(0)||ai||i a ui 

-A-(^(t))||o a ||i a «i, 
fo 22 = id t u 2 + -Au 2 - i (jte, V 2 V (i, q- 2 (t)) 2/2) w 2 - «"(0)||o2||i a u2 

-A-(-«g(i))||oi||i a U2. 

The last term in each expression accounts for a coupling, revealing a leading order inter- 
action of the two wave packets. This coupling can be understood rather explicitly, since 
it consists of a purely time dependent potential. Solving the equations b 2 j = 0, we infer, 
with obvious notations adapted from ( 11.5b , 

2 



ui(t,i/i) l/i)«q> ( -iiif(0)||ai||i a -*||a 2 |||, 2 J K(Sq(s))ds 

u 2 (t,y 2 ) =u 2 in (t,y 2 )exp(-itK(0)\\a 2 \\ 2 L2 -i\\a 1 \\ 2 L2 f K{-Sq(s))ds 



The presence of these phase shifts accounts for nonlinear effects at leading order in the 
approximate wave packet tp^ pp : nonlinear effects in the case of a single wave packet, and 
nonlinear coupling, since we assume ||aj||x,2 7^ 0. For the remainder terms, we have the 
(rough) pointwise estimate 



(2.7) \r*{t, yj )\ < C{T)J-e{ yj f\ Uj {t, yj )\ ( 1+ £ |K(t)||| 2 ] , te [0,T 

fe=l,2 



The remainder is the sum of the term s r e ^ v and £ a+3 ^ 2 (^fc * \ u k\ 2 )uj, k = 1,2, so 
this estimate is an easy consequence of ( 12.3b and (12.6b . To be precise, this estimate is 
valid up to the rectangle terms that we have discarded so far, when we have developed 
(K * |^appj 2 _)' ( /app- ^ e will see in Section [3] that they satisfy a similar estimate (see 
Corollary 13. 2b . 

2.3. The case a = 1/2. We still have b 0j = b^f, but now 

hi = -i (gi(t) -Pi(t)) ■ Vui - yi ■ (pi(i) + W (*, qi {t))) m 

-X(0)|| ai || 2 2Ul -if (<5g)||a 2 ||| 2 ui, 
612 = -i feW - ]*(*)) • Vu 2 - y 2 ■ (pj,(i) + W (t, $,(«))) "2 

- (0)||a 2 ||i 3 u 2 - ^ ||oi||i a « 2j 
631 = tftui + -Atti - - (yi, V 2 V{t, qi (t)) yi) ui - ||ai||i 2 yi ■ VK(0) Ul 



\a2\Wm ■ VK(dq) Ul + Vtf (0) • Gi(<)mi + VK(Sq) ■ G 2 {t)u x , 
u 2 + ~Au 2 - ~ (iftj, V V (i, ga(t)) 2/2) «2 - \\a 2 \\ 2 L 2y 2 ■ VK(0)u 2 
\a 1 \\l^-VK{-5q)u 2 + VK{Q)-G 2 {t)n 2 + VK{-8q)-G 1 {t)u 2 . 
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The first line in bij is zero if (qj,Pj) is the classical trajectory given by (II. 3t . On the 
other hand, it does not seem to be possible to cancel out the second line in bij, even by 
modifying ( 11.31 i: we have three sets of terms, involving Viti, y\U\ and u\, respectively, so 
they must be treated separately. As in (£), we then modify the general strategy, and allow 
boj to depend on e. We alter the hierarchy as follows: 

6gi - -«i (Si (t) - Pi (*) • qi (t) + + V (t, 9i (t)) 



"02 
&12 



+ ^(0)|| fll || 2 L2 + V^K(Sq) \\a 2 f L2 j, 
-U 2 (s 2 (t) - p 2 (t) ■ q 2 (t) + + V (t, q 2 (t)) 



+ ^K(0)\\a 2 \\ 2 L2 + JiK (Sq) ||ai||| 2 J , 
-i ( qi (t) - Pl (t)) • Vu x - 2/i ■ (j>i{t) + W (t, qi (t))) ui, 
-i (q 2 (t) - p 2 {t)) • Vu 2 - y 2 • (pa(t) + W (t, ©(*))) u 2 , 



and we leave unchanged. Like before, bij = provided that (qj,Pj) solves ( 11.3b . The 
novelty is that we now consider modified, e-dependent, actions: 



sm = 



\pi(s)\ 2 -V(s ) qi(s)))ds 



-tv/iir(0)||ai||| 2(R<i) - v^lKHi^Rd) / K(8q(s))ds, 

J 



s 2 e W = 



-| P2 ( s )| 2 -y( S7(Z2 ( s )))d s 



- tVeK(0)\\a 2 \\ 2 L 2 {Rd) - v / e||ai||i2(Rd) J if (Sq(s)) ds. 

These expressions are exactly those given in the introduction (|1.9l l. The equations b 2 j = 
are envelope equations, which are nonlinear since Gk is a nonlinear function of Uk- Note 
however that Gk yields a purely time-dependent potential. Consider the solution to 

id t u x + ^Aui = i (yi, \7 2 V (t, qi[t)) Vl ) Ui + \\aif L2 yi ■ VK(0)ui 

+ \\a2Wlm ■ VK(Sq)ui, 
id t u 2 + ^Au 2 = 1 (y 2 , \7 2 V (t, q 2 {t)) y 2 ) u 2 + \\a 2 \\ 2 L2 y 2 ■ VK(0)u 2 



+ ||oi||i>ife- VJC(-5g)u a . 



Set 



iti(i, yi) = Gi(t,yi)exp i 



u 2 (t,y 2 ) = u 2 (t,y 2 )exp 
where 



(vJf(Q) • Gi(s) + Vif • G 2 (a)) *) , 

(V-fT(O) • G 2 (s) + VK (-5q(s)) • j . 



R<< 
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It is clear that \v,j\ = hence Gj = Gj, and so u\ and u 2 are such that 621 = &22 = 0, 
and correspond to the envelopes introduced in £11.41 Finally, we still have a remainder term 
satisfying ( 12.71 i (up to the terms treated in 

2.4. The case a = 0. Now all the coefficients bij are affected by the nonlinearity: 

601 = -u x (s x {t) - Pl (t) ■ &(t) + + v («,«!(<)) + ^(o)|K||| 2 

+ IK III*), 

602 = -u 2 (s 2 (t) - P2 (t) ■ q 2 (t) + + V(t,<fc(t)) + ^(0)1108 Hi, 

2 A 

\L 2 h 



+ K{~5q)\\ ai \ 2 

b u = -i (qi(t) - pi(t)) • V«i - 2/1 • (pi(t) + W (t, ?i(t))) «i - HoiH^yi • VK(0)ui 

- IMI^yi • Vif Ul + VJf(O) • G?i(t)«i + VK (8q) • G 2 (i)«i, 

&12 = -*(&(*) -P2W) ■ Vu 2 -2/ 2 • (R(t) + W(t,? 2 (t)))«2 - lloalliaKa ■ VK(0)u 2 

- Hoilliajft • Vif (-tfg) u 2 + Vtf (0) • G 2 (i)w 2 + Vif (Sq) ■ G x {t)u 2 , 

hi = id t U! + ^A Ul - i (yi,Mi(i)i/i) «i + (V 2 ^(0)Gi(t), yi) ui 

+ (V 2 /^)G 2 (t), j/i > (z, V 2 K{0)z) |«i (t, z)| 2 daA mi 

& 22 = i&tta + ^Au 2 - ^ (y2,M 2 (t)y 2 ) u 2 + (V 2 K(0)G 2 (t),y 2 ) u 2 

+ (TpKi-SqlGi®, y 2 ) ^~\[J ( 2 > V 2 K(0)z) \u 2 (t, z)\ 2 dz^j u 2 

where we have denoted 

Mi(t) = ||ai|li 2(R£i) V 2 ^(0) + \\a 2 \\ 2 LH ^ 2 K (6q(t)) + V 2 X V (t, qi (t)) , 
M 2 {t) = \\a 2 \\ 2 L2{Tld) V 2 K{0) + ||ai||i 2(Rd) V 2 ^ (Sq(t)) + V 2 X V (t, q 2 (t)) . 

Similar to the case a = 1/2, we incorporate the last term of bij into boj, that is we modify 
the action as follows: 

Sf(t) = J* (\\pi(s)\ 2 - V (s, qi (s)) - K(0)\\"i\\h K (*«(«)) 
+ y/eVK(0) • Gi(s) + VeVlf (5q{s)) ■ G 2 (s)W 

Slit) = J* (\\P2(s)\ 2 - V (s, q 2 (s)) - K(0)\\a 2 \\ 2 L 2 - K (-Sq(s)) |K | 
+ VeVA"(0) • G 2 (s) + v^VJsT • Gi(s)W 



2 

L 2 
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Note that for 5| to be well defined, we have to first define Uj, for which we solve the 
envelope equations, given by b 2 i = ^22 = 0. Canceling the terms b\j yields the modified 
system of trajectories: 

( qi{t)=pi(t), 

Mt) = -VV(t, qi (t)) - || ai || 2 L2 VX(0) - \\a 2 \\ 2 L2 VK ( qi (t) - q 2 (t)) , 

92(f) =P2(t), 

{ Mt) = -VV(t,q 2 (t)) - ||a 2 ||| 2 VX(0) - llailliaVA- (q 2 (t) - gi (t)) , 

which is exactly dl . 12b . The remainder term still satisfies (12.71 > (up to the terms treated in 
SJ3]). We will examine more carefully the envelope system in Sj5] 

3. Analysis of the rectangle interaction term 

In the previous section, we have left out the rectangle terms, claiming that they are 
negligible in the limit s — s- 0. In this section, we justify precisely this statement. Since the 
two terms that we have discarded are similar, we shall simply consider the first one: 



2e- d / 2 Re J K(z)e«*-+*X t >*-*V e u 1 [t, Vl - -Z= , { ,. : „, , _ _ 

Notice that we have not yet expressed the phases <fik in terms of the variable yi, and that 
the expression of (f>k varies according to a — 1, a = 1/2, or a = 0. We shall retain only a 
common feature though, that is, (f> e k {t, x) — Of (i) + x ■ Pk(t), where only the purely time 
dependent function may depend on e (when a 6 {1/2, 0}), and the spatial oscillations 
are singled out. Since x = qi (t) + \fey\, we get, once the real part and the time oscillations 
are omitted: 



u 2 t, yi + 



6q(t) 



dz. 



dz. 



Changing the integration variable, and introducing more general notations, we examine: 



(3.1) 



r(t, yi )= JlC - z)) e^'W/^ti! (t,z)u 2 (t,z+ 6 -^) 



dz. 



The main result of this section is stated as follows. 

Proposition 3.1. Let T > 0. Suppose that JC £ W*>°°, u G C([0, T]; E fe ) with fe,leN, 
and consider I e defined by d3.lt . There exists C > independent of e € (0, 1], /C, Mi and 
u 2 such that 

sup ||,r e (t, -)lk^ (R . d) =s; c* e mi "^^>/ 2 ||/c|| M ^^||^ 1 || £oo([0 , r];Sfc) ||^ 2 || iac([0 , r]:Sfc) . 

te[o,T] 

In view of the computations performed in Section [2] this result has the following con- 
sequence. 

Corollary 3.2. Consider ip £ pp given by ( 12.11 ), derived in Section]2\ whose exact expression 
varies according to the cases a = 1, a = 1/2 or a = 0. Let T > 0, and suppose 
uj e C([0, T]; S 3 ). Then ^ pp € C([0, T]; £ 3 ) satisfies ip £ a 



ipf t=0 and 



'app ~ U"J J) " J ^app|t=0 
feft^pp + Y A^ pp = V (t, l) V a pp + >/i (A" * IV'appI 2 ) ^app + ^ 

where r £ £ C(R+; L 2 (R d )) is smc/z f/iaf there exists C independent of e with 

sup IKlXdj^Rd) < CVe- 
te[o,T] 



18 



R. CARLES 



Remark 3.3. At this stage, the property Uj G C([0,T];£ 3 ) is established in the cases 
a = 1 and a = 1/2. It will require some work to prove it in the case a = 0; see Section|5] 
The assumptions of Corollary 13. 2l are fulfilled, modulo the proof of Proposition[ 



Remark 3.4. Proposition B.ll is a refinement of IU21 Proposition 6.3], in the sense that the 
power of £ on the right hand side is as large as we wish, provided that JC is sufficiently 
smooth, and that the functions u% and u 2 are sufficiently localized in space and frequency. 
Identifying precisely the norms of JC, u\ and u 2 , involved in order to get such an error 
estimate, will turn out to be crucial to prove Theorem 11.91 at the level of the bootstrap 
argument presented in Section|7] 

3.1. A microlocal property. The proof of Proposition B.ll is based on the following re- 
mark: the function that we integrate is localized away from the origin in phase space: 

Lemma 3.5. Suppose (qio,Pio) 7^ (<?20jP2o)- I n either of the cases a = 1, a = 1/2 or 
a = 0, the following holds. For any T > 0, there exists rj > such that for all t 6 [0, T], 

\Sq(t)\^T], or\Sp{t)\^ V . 

Proof. We argue by contradiction: if the result were not true, we could find a sequence 
f„ G [0,T] so that 

\5q(t n )\ + \6p(t n )\ — ► 0. 

n— >-oo 

By compactness of [0, T] and continuity of (qj,pj), there would exist € [0, T] such that 

qi(t*) = q2(U), Pi(U) = P2(t*)- 

In the cases a = 1 and a = 1/2, (qj,Pj) is given by the classical Hamiltonian flow ( 11.31 ): 
uniqueness for ( 11.3b implies (<Zio,Pio) = (<720)£>2o)> hence a contradiction. 

The case a = is a bit more delicate, since (qj , pj ) is no longer given by a Hamiltonian 
flow. From ( 11.121 1, we infer: 

f d(6q) 



dt 
d(Sp) 
dt 



Sp, 



- W (t, q 2 (t)) - W (t, Ql (t)) + \\ ai \\ 2 L2 (VK (Sq(t)) 



\a 2 f L 2 (VK(0)-VK (Sq(t))) 



In view of Assumption ll.il there exists C independent of t such that 

^C(\Sp\ + \Sq\). 



d(Sq) 


+ 


d(S P ) 


dt 




dt 



Gronwall's Lemma yields a contradiction, and the lemma is proved in the three cases. □ 

3.2. Proof of Proposition 13.11 From Lemma [331 if suffices to prove the estimate of 
Proposition l3.1l in either of the two cases |<5g(t)| ^ 77, or |<$p(t)| ^ rj. 
First case. If |<Sg(t)| ^ rj, we use Cauchy-Schwarz inequality to infer 



\F(t,y)\^\\fC\\ L , 



^ \\IC\ 



"I K(M)I 



z + 



Sq(t) 



Z + 



M0 



«2 



dz 



\ui(t)\\ S k\\u 2 (t)\\ S h sup (z)~ 
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In view of Peetre inequality (see e.g. ll2ll4Tl'). 

Second case. If \Sp(t) | ^ r\, we perform repeated integrations by parts (like in the standard 
proof of the nonstationary phase lemma, see e.g. 0) relying on the relation 

Note that since we assume JC G W '°° and Uj G £ fe , we perform no more than min(£, k) 
integrations by parts, and Cauchy-Schwarz inequality yields 

\F(t, y )\ < ^\m wt ,o \\u l m&\\Mm&^ Mt,k)/,> - 

The proof of the proposition is complete. 



4. Proof of convergence in the critical case 



In this section, we complete the proof of Theorem 1 1.3 1 First, we recall that as a conse- 
quence of 1171 fill , the system for the envelopes in the linear case is well-posed in E fe : 

Lemma 4.1. Let k G N, and a G Yi k . Then ( 11.51 ) has a unique solution u G C(R+; 
In addition, the following conservation holds: 

l|w(*)IU a (R<*) = ||a|| L 2 (R£i) , Vi ^ 0. 

We infer that if a\, a 2 G S 3 , then u 1: u 2 , given by (11.81 . belong to C(R + ; £ 3 ). Corol- 
lary [3T2]iiiiplies that ipl pp satisfies 

e 2 

ied t *p s app + = ^app + e{K* |^ pp | 2 ) ^ pp + er £ ; V apP (0, x) = V> £ (0, as), 

where the source term r £ satisfies: 

VT>0, 3C7>0, sup \\r E (t)\\ L 2 {Rd) ^ Cy/e. 
fe[o,T] 

Denote by w £ — ip £ — ipl pp the error term. It satisfies 

tsd t w £ + = Vw £ + e ((K * \^\ 2 ) ^-{K* |^ pp | 2 ) ^ pp ) - sr £ , 

with it>f t _ = 0. Writing 

(if * |V> £ | 2 ) 1p £ - (K * |^app| 2 ) V'app 

= (if * \W £ + ^ pp | 2 ) (lif + ^ pp ) - {K * |^ pp | 2 ) </>app 

= (if * | W £ + ^ pp | 2 ) W £ + (K * (\W £ + Vappl 2 " IV'appI 2 )) V>app, 
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energy estimates yield, for t € [0, T], and since ip £ , tp^pp (hence w £ ) are uniformly bounded 
inL°°(R + ;i 2 (R d )): 

\\vf{t)\\» < T * (K + ^appl 2 - l^appH W|| LOO ll^ppWIlL-dfl 

Jo 

+ / \\r £ ( S )\\ L 2ds 
Jo 




for C > independent of e € (0, 1] and t ^ 0. Theorem 1 1.3 1 is then a consequence of 
Gronwall's Lemma. 

Remark 4.2. Assuming that we have proved the property Uj € C([0, T]; E 3 ) in the case 
a = 0, which will stem from Proposition |1.81 the conclusion of Corollary 13. 2| holds. How- 
ever, the estimate given by the above approach is not satisfactory in the cases a = 1/2 and 
a = 0. We could prove this way: 

W £ (t)-ip s app (t)h> <<Vie c */ £l ~ a , te[0,T], 

for a = 1/2 and a = 0, respectively. Contrary to the case a = 1 (where Gronwall's 
Lemma yields a similar estimate), we can only conclude that ip £ — ipl pp is goes to zero on 
a small time interval: there exist c > and 8 > such that 

sup ||^W-^pp(t)IU=— >0. 

O^Kce 1 -"! lne| 8 

Corollary 13. 2| is a consistency result, which is not enough to infer convergence. This can 
be understood as a feature of supercritical regimes: a different approach is needed, which 
requires more regularity from V, K, and the initial data cij. 



5. The envelope equations in the case a = 



(5.1) 



In this section, we prove Proposition 11.81 We first remark that the last two terms in 
each equation involved in ( 11.13b correspond to purely time-dependent potentials, and can 
be treated thanks to the gauge transforms 



«i(t, y) = Ui(t, y) exp ^ — i j j (z, V 2 K(0)z^ \u±(s, z)\ 2 dzds 

(z, V 2 K (Sq(s)) z) \u 2 (s, z)\ 2 dzds^ , 
(z,V 2 K(0)z) \u 2 {s,z)\ 2 dzds 
- / I I (z,V 2 K (-5q(s)) z) \ui(s, z)\ 2 dzds). 



V2(t,y) = u 2 {t, y)cxp ( 
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Since K is real-valued, we have \vj(t, y)\ = \uj(t, y)\, and (15.11) is equivalent to 

l (t,y)=v 1 (t,y)eiq)(i ^ f (z 1 V 2 K{0)z)\v 1 (s 1 z)\ 2 dzds 

(z, V 2 K (Sq(s)) z) \v 2 (s, z)\ 2 dzds) , 



(5.2) 



u 2 (t,y) = v 2 (t,y)exp(i j j (z,V 2 K(0)z) \v 2 (s, z)\ 2 dzds 

(z, V 2 K (-5q(s)) z) \vi(s, z)\ 2 dzds 



+ i 



Formally, (ui, u 2 ) solves (11.131 > if and only if (vi,v 2 ) solves 

\(V,Mi(t)y) Vl -(V 2 K 
(\7 2 K(Sq(t))G 2 (t),y) Vl , 
[ -(y,M 2 (t)y) v 2 -(W 2 K(0] 
(W 2 K(-Sq(t))G 1 (t),y)v 2 , 



(5.3) 



id t vi + ^A Vl = i (i/,Mi(t)i/> vx - (V 2 K(0)Gi(t),y) Vl 



id t v 2 + X -Av 2 = i (y,M 2 (t)y) v 2 - (V 2 K{Q)G 2 {t), y) v 2 



with the same initial data, Vj\ t -Q = a,j, j = 1,2, where the bounded, symmetric matrices 
Mi and M 2 are defined in (11.15b and (I1.161 l. respectively, and where we have kept the 
notation 

Gj(t) = J z\ Vj (t,z)\ 2 dz. 

Note that the terms involved in the gauge transforms are well defined when the functions 
are in E fc with k 1, so Proposition 1 1 . 8 1 stems from the following: 

Proposition 5.1. Let (qi,pi, q 2 ,p 2 ) be given by Lemma \L6\ and a±,a 2 G S fc with k 1. 
r/ien ( 15.3b /las a unique solution (vi,v 2 ) G C(R+;S fc ) wif/z initial data (01,02). /« 
addition, the following conservations hold: 



M*)ll 



L 2 (R d ) 



£ 2 (R d), Vt ^ 0, j = 1,2. 



Proof. The main difficulty is that since the last two terms in each equation involve time 
dependent potentials which are unbounded in y, they cannot be treated by perturbative 
arguments. So to construct a local solution, we modify the standard Picard iterative scheme 
in the same fashion as in [9], to consider 

' id t v[ n) + \Av[ n) = \ (y,M 1 (t)y) - (v 2 if (0)Gf "^(t), j,) ^ 

-(v 2 K(Sq(t))G ( r 1} (t),y)v^\ 

id4 n) + ±Av™ = \ (y,M 2 (t)y) »W - (v 2 K (0)^ (t), y) «£° 

-(v 2 tf (-*?(<)) Gf-^y)- 



(5.4) 



An) 



with Vju—Q = a j f° r a H n 7 v f\ti v) — a j(y)' ar, d 



Gf\t) 



1 * 


( fc ) 1 \ 


2 

ds 
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At each step, we solve a decoupled system of linear equation, with time dependent poten 
tials which are at most quadratic in space. If G± , G 2 n € LJ^ C (R+), ifTHl ensure: 
the existence of € C(R+; L 2 (FL d )). In addition, we have 



L 2 (R d ) 



| L 2 (R(i ) Vi > 0, j = 1,2. 



Applying the operators y and V to ( 15.4b yields a closed system of estimates, from which 
we infer that e C(R+; S), hence Gp € Lj^ c (R+). Therefore, the scheme is well- 
defined. Higher order regularity can be proven similarly: for k ^ 1, by applying k times 
the operators y and V to (15. 4K we check that i>J"' ) , w^™'' € C(R + ; E fe ). As a matter of fact, 
due to the particular structure of (15.4) . the only informations needed to prove this property 
are aj € S fc and vf~ 1] e C(R+; S). 

To prove the convergence of this scheme we need more precise (uniform in n) estimates. 
A general computation shows that if v solves 

idtv + i Av = - (y, M (t)y)v + F(t)-yv, 

where M(t) is a real-valued, symmetric matrix, and F(t) is a real-valued vector, then 
G(t) — / z\v(t, z)\ 2 dz satisfies formally 



G{t) = Im / vVv =; J(t), 



J(t) 



(M(t)y + F{t)) \v(t, y)\ 2 dy = -M{t)G{t) - F(t)\\v\ 



L 2 ' 



where the last expression uses implicitly the fact that the L 2 -norm of v is independent of 
time. We have in particular: 

G(t)+M(t)G(t) = -\\v\\l 2 F(t). 



In our case, this yields: 

(5.5) G < ? ) +M 1 {t)Gf ) = ihllz^/vdnr;', 



(5.6) G 
Let 

We have 



{n - 1) + \\a 2 \\l 2 W 2 K(6 q (t))G ( r 1 \ 
(") i us t+\ri( n ) — iu I|2 yv2 isKwnt™- 1 ) i n„ i|2 t-iIts < a^i+w n^ n ~ x ) 



M 2 {t)G\ 



f n (t)= G[ n \t) + G 2 n '(t) + G[ n >(t) + G\ n >{t) 



\a 2 \\i 2 V'K(0)G 2 



: (n), 



\ ai \\i^K(-5q(t))GT 



f n {t) < 2 £ (|G(. n) w| \Gf\t)\ + |c< n) (f)| \G { ;\t) 

3=1,2 



for some C independent of t and n, since V 2 V, V 2 K 6 and where we have used 
d5 .5b -( T5~6b and Young's inequality. By Gronwall's Lemma, we infer 

fn(t) < f n (0)e ct + C f e c( *- s )/ n -i(s)ds. 
Jo 
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With our definition of the scheme, /„ (0) does not depend on n: 



fn(0) = E ( Im / 
.7=1,2 V J 



z\a,j{z)\ 2 dz 



= :C . 



Therefore, 

/»(*) ^ C e ct + C [ e c ^ f n -i(s)ds, 
Jo 

and by induction, we infer 

f n {t) < 2G e 3Ct , t > 0. 

By using energy estimates (applying the operators y and V successively to the equations), 
we infer that there exists C\ independent of t ^ and n such that 



J=l,2 



< Cie 



Cit 



The convergence of the sequence , «2 ) then follows: we check that v n converges 
in C([0,T]; E) if T > is sufficiently small. To simplify the presentation, we present 
this argument in the case of a single envelope equation, the case of J5.4I ) bearing no extra 
difficulty: 

(5.7) idtv^ + 1a»W = \ (y,M(t)y)vW + (qWG*"-^*), !/)«<">, 
where Q(t) is a real-valued, symmetric matrix, with Q E L°°(R + ). Denoting by 

H(t) = -±A+~(y,M(t)y), 



we have 

i 



+ (Q(t) (g^^-G^ 2 ^)) ,y)« (n " 1} 
Energy estimates and the above uniform bound yield 



v (n)^_ v (n-l)^ 



L' 1 



<G 



G ( "- 1) (s)-G ( "^ 2) (s) 



f/.s 



^ G 



By applying the operators y and V to ( 15.71 ), we obtain similarly: 



tWfi)-^- 1 '^) 



< G 



v ( ™- 1) (s)- U ^-^(s) 



(n-2) , 



ds. 



Therefore, we can find T > such that the sequence i;W converges in C([0,T]; £), to 
u G G([0,T]; The uniform bounds for the sequence v^ n ' imply that v is global in 
time: v £ G(R+; S fe ), with E fc -norms growing at most exponentially in time. □ 

6. Convergence in supercritical cases: scheme of the proof 

We present the proof of Theorem 1 1.9 1 in details; the proof of Theorem l 1.5 l ean easily be 
adapted (see Remark Rx2l below). 
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6. 1 . The general picture. In 1 3 9 1, where the case of only one wave packet is considered, 
the proof of stability relies on a change of unknown function: writing 

r(t,x) = s - d / i u e U^A 

with S £ , q and p as given by the construction of the approximate solution, it is equivalent 
to work on ip £ or u £ in order to prove an error estimate, since 

Wit) - ^pp(*)IU»(R*) = - «(*)IU'(R<)- 

Passing from the unknown tp £ to u E amounts to using very fine geometric properties related 
to the dynamics: the modified action S £ , and (q, p). One changes the origin in phase space, 
to work in the moving frame associated to the wave packet. In the case of two wave packets, 
there are two moving frames, so the approach that we follow is different. We construct a 
solution to ( II . lb of the form 

(6.1) = £- d/i T u'Jt, X ~ q / i {t) ) e <sj(t)+p 3 .( t ).(*- 9i ( t )))/e 5 

where the quantities (qj,Pj) and S £ are those given by the construction of tp^ pp , so we 
consider two unknown functions, u\ and To do so, we derive formally a system for 
(u\ , which is morally equivalent to ( ll.lt : rigorously, the solution to this system yields 
a solution to dl.lt , and by uniqueness for ( II. j} , the relation ( 16. It is valid. In turn, the 
construction of the solution (u\ , u?,) on arbitrary time intervals relies on a bootstrap argu- 
ment, consisting of a comparison of a modification of {u\ , itf) with (ui,W2), defined in 
(11.13b . This modification eventually corresponds to the presence of the phase shifts Oj in 
Theorem ll.9l 

In order to shorten the formulas, we consider indices in Z/2Z: typically, qj stands for 
<7i whenever j = 1 or 3. Plugging ( 16. U into ( 11.11 ) in the case a = 0, we find: 

F 2 

ied t tp £ + —Aif; £ - Vtp £ - (K * \*fj £ \ 2 ) iP £ = e- d/i e l ^ {t ' x) N £ , 

where we have denoted 

<P £ j {t,x) = S](t)+p j (t)-(x-q j {t)), 

N £ = isd t u £ - u £ d t <P £ + -Auf - ^^«| - K (t, q 3 {t) + y jy /S) u £ - Vf L ^, 
Vf L (i, % ) = J K (Vife - z)) \u%t,z)\ 2 dz 

+ J K (<U - Qj+i + VsiVj - z)) \u e j+1 {t, z)\ 2 dz 

+ 2Ree i ( SE j~ s J+ 1 ~ qj ' Pj+qj+1 ' Pj+1 ^ £ x 

x J K - z)) e iz <^- p ^/^u%t, z)u £ +l (t, z + qj ~j! +1 ^ dz. 

Note that the computations which we do not detail correspond to the computations pre- 
sented in Section|2l up to the fact that now, we do not perform Taylor expansions for V or 
K. As in Section|2l we distinguish the variables y\ and i/2- 

Our approach consists in considering the set of coupled, nonlinear equations 



\ e i(S e (t)+(x-g(t))-p(i))/£ 
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It is important to notice at this stage of the construction that this system conserves formally 
the L 2 norms: since we naturally impose Uj* t _ = a 3 , we have 

ll u i(*)llL 2 (R d ) = l|a?IU 2 (R< i )> 3 = 1^2, 

as long as (u\ , m|) is well defined. This property is the reason why we can perform impor- 
tant reductions in the system. Taking into account the expression of the modified actions 
Sj, we find: 

N] = ied t u) - Vjj (s?(t) + y/ipj(t) ■ y 3 - Pj (t) ■ q 3 {t) 



: IE 



^^u] - V (t, q 3 (t) + y iy /e) u) - Vf L (i, Vi )u) 



d t u) - (Vepj{t) ■ y 3 - pj(t) ■ qj (t)) u) 
- (\\Pj\ 2 - V(t, qj ) - K(0)\\ aj \\ L2 - K(q 3 - qj+1 )\\a j+1 f L ^ u) 
+ (V^Vif(O) • Gj{t) + ^feVK{q 3 - q j+1 ) ■ G 3+1 (t)) u) 



- V (t, q 3 (t) + y jy /e) u) - V 3 NL (t, Vj )u) 



= ied t u 3 + -Au^ - Vep 3 (t) ■ y 3 u) 
-(V(t, qj (t) + yj V^ -V(t,2,-))«S 

- (v^(t, Vj ) - K(0)\\a 3 \\ L2 - K(q 3 - <z, + i)K+i||i 2 ) u| 
+ (ViVX(O) • G,(t) + v/iV^fe - qj+1 ) ■ G 3+1 {t)) u). 
If we now take into account the expression of p 3 , given in (I1.12l i. we infer: 

NI=isd t ue+ e -& U * 

+ sTeyj ■ (yV{t, qj ) + \\ aj f L2 VK(0) + \\a 3+1 \\ 2 L2 VK( q3 - q 3+1 )) 
-(V (t,q 3 (t)+y 3 ^) -V(t, qj )) u) 

- (v 3 Nh (t,y 3 ) - K(0)\\a 3 \\ 2 L2 - K(q 3 - %+ i)||a J+1 || 2 L2 ) u] 
+ (y/2VK(p) • G 3 (t) + V^VK( qj - q 3+1 ) ■ G 3+1 (t)) u*. 
It is now natural to introduce the following notations: 

Vf(t, y 3 ) = ~ (V (t, q 3 (t) + y jy /i) - V (t, q 3 (t)) - Jiy 3 ■ W (t, q 3 (t))) 



K 



j,diag 



(t, y 3 ) = - (K (Vijy) - K(0) - ^ey 3 ■ Vif(O)) , 



Kj,off(ti Vi) = -( K («?' _ *+i + y/eVj) - K(q 3 - q 3+1 ) - ^fey 3 ■ VK{q 3 - q 3+1 )) . 
From the assumptions on V and K, there exists C > independent of e £ (0, 1] such that 



(6.2) 



E ||W*) 



E 

2<|a|<6 
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where Kj stands for di or Kj oS , indistinctly. In view of Taylor's formula, we have: 

(6.3) Vf(t, Vj ) = f ( Vj , V 2 V (t, qj(t) + Byjy/e) y s ) (1 - O)d0, 

Jo 

(6-4) #J, diag (*, Vi) = I (y 3 ,V 2 K (e yj V£) Vi) (1 
(6.5) 

Therefore, we consider the coupled system (coupling is present through Kf): 



K e j,os(t, Vj) = / (Vi, V 2 K ( qj (t) - q j+1 (t) + 9y 3 V^) Vj) (1 



(6.6) { 



id t u) + -An* = Vf(t, + (Kl diag * \uf) u| + (K] <oS * |u* +1 | 2 ) 
~VK(fl)- (J z(\u^z)\ 2 ~\u,j{t,z)\ 2 ) dz^j 
-^fK{q 3 -q 3 + 1 )- (J z(\u £ J + 1 (t,z)\ 2 -\u J+1 (t,z)\ 2 )dz^j 



(2ReW*(t, yj )) u), 



j(S^-S^ +1 -q j -p j +q j + 1 -p j + 1 )/e 



x J K (V~e(y 3 - z)) e^-^y^it, z)u^ +1 (t, z + gj 



dz. 



with 



6.2. Further simplification and bootstrap argument. The last three terms in (16.6) are 

singular in the limit e — > 0. However, the singularity of the last term is expected to 
be artificial, since in view of Proposition 13.11 it should even be small, provided we have 
uniform estimates for Uj in S 3 . The other two singular terms have an interesting feature: 
they are real-valued, and depend only on time, so we can treat them thanks to a gauge 
transform. Introduce 



1 



m s j + 5 Au| = Vffryj)^ + (iq diag * m 2 ) v? s + (Kl oS * |^ + il 2 ) «? 



(6.7) 



with initial data u e ,. „ = a 7 , and where we have denoted 



;( 



2ReWf )m, 



W* = e 



l ( S j- S 'j + l-^-W+?3 + l-W + 1 )/ £ e i ( e 3- e i + l)- 



x / K (Ve(» 3 - - z)) e **-Cw-w+OA^( t) z) u] +1 [ t)Z + 



dz, 



with 



0f(f) = ^=J\k(0) ■ (J z{\u]{s,z)\ 2 - \u 3 (s,z)\ 2 ) dz^j ds 

We then have: Uj(t, y) = Uj(t, y)e t6 ^ t \ Note that = \u^\, so it is equivalent to pass 
from Uj to Up or from u £ j to u|. In view of these reductions, in a first approximation, 
Theoreml 1 .91 stems from: 
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Theorem 6.1. Let d ^ 1 and a\,a2 G E 6 . Assume that V and K are real-valued and: 
V E C 6 (R+ x R d ; R), and d?V £ L°° (R+ x R d ) , 2 < < 6. 
^ e vK 6 '°°(R d ;R). 

Lef T > 0. There exists Eq > smc/z that for e £ (0,£o], ( 16.71 /icw a unique solution 
(u\ , m|) S C([0,T]; E 3 ) 2 . Moreover, there exists C independent of e G (0, £o] 5Mc/z f/iflf 

(6.8) sup ||«!(t)-ui(t)||j33+ sup ||fii(t)-ua(*)|| E s <CVe. 

te[o,T] te[o,T] 

Several comments are in order. First, this result implies that for j = 1,2, is bounded 
on [0, T], uniformly in e. To get the result stated in Theorem 1 1.91 we will prove that the 
functions 0| converge as e — > 0, by performing a second order asymptotic expansion of 
(Theorem l6. 1 I vields the first order asymptotic expansion). 

Even in the case of a single wave packet, this result is new, since we do not assume 
VJC(O) = 0. In that case, the last term in (16. 7t vanishes, and the proof that we present 
below becomes simpler. 

The proof is based on a bootstrap argument detailed in Section[7] For fixed e > 0, d6.7t 
has a unique, local solution: (ufjii!) G C([0, t e ]; E 3 ) 2 , for some t £ > 0. This can be 
proven by adapting the approach presented in Section To prove the theorem, we use 
energy estimates to prove that so long as (u\ , {t|) is bounded in C([0, r]; S 3 ) 2 , t £ ^ r ^ 
T, j6T8l ) is true. Therefore, choosing e > sufficiently small, (wf,<i/f) € C([0,T]; E 3 ) 2 
for e G (0, So], and d6.8t is satisfied. 

The reason why we work in E 3 and not in a larger space is that we want to be able to 
neglect WJ: because of the singular factor 1/e in front of the last term in ( 16.7) , we need to 
prove W? = o(e), and Proposition |3.1| suggests that we need to work in E 3 , in which case 

= C(e 3/2 ). To differentiate Vf and K £ three times (we work in E 3 ), d63T>— (|63T> and 
Proposition 13 . 1 1 suggest to work with the regularity stated in Theorem |6T| (the same as in 
Theorem |1.91 >. 

Remark 6.2. In the case a = 1/2, one can consider that all the terms involving K are 
multiplied by yfe. As a first consequence, it is enough to work in E 2 to prove that the 
term involving Wj is negligible. By working in E 2 , we only need to differentiate and 
Kj twice, hence the regularity assumption in Theorem 1 1.5 1 Finally, since the phase shift 
relating u 6 and u £ is multiplied by y/e, it is 0(y/e), as opposed to 0(1) in the case a = 0. 

7. The bootstrap argument 

In this section, we prove Theorem |6.1| More precisely, we focus on J6.8t , in view of the 
discussion at the end of Section [6] 

In Section [6] we have essentially resumed the computations of Section [21 U P to two 
aspects: 

• We have not used Taylor's formula for V and K. 

• The terms 9 e do not appear in the case of the Uj's (replacing u e with Uj in the 
expression of 9 e yields 9 e = 0). 

In Section [3] we have seen that the analogue of the term Wj (or, equivalently, WJ) is 
negligible in the limit e — > 0. These properties can be summarized as follows: the functions 
iti and U2 solve 

(7.1) id tUj + -AUj - VfUj + (^J.diag * Wj I') U, + (K? off * \U J + 1 1 2 ) Uj + p), 
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where p e is given by the formula: 

P) = (K - Vf) Uj + ((Kl dias - if° diag ) * \ Uj \ 2 ) Uj 
+ ((Kl oS -Kl oS )*\u J+1 \ 2 )u J , 

where V®, K® are given by (16.3b — (f6T5b with e = 0. We infer from (16.31 — (16.5b and Propo- 
sition l 1 .8l that for all T > 0, there exists C > independent of e 6 (0, 1] such that 

sup 11^0)11 S 3 C^fe. 
te[o,T] 

Since the bootstrap argument runs in £ 3 , it is natural to work with such an estimate for 
the source term. This in turn imposes to work with a,j E E 6 , as well as V and K as in 
Theorem|6.1| 



Set w e - 



f subtracting d7.lt from (16.7) , we see that the error satisfies the coupled 



system, for j = 1,2, 



(7.2) 



1 



idtw? + -Awl = VfvJj + (if? dlag * |^ £ | 2 ) u) - (if|, dlag * hf ) Uj 

+ {Kl oS * |^ + i| 2 ) u e j ~ (i^off * K'+i I 2 ) «j 
1 



2 Re W? Hi 



Pi 



with initial data w^, t _ Q = 0. Fix T > once and for all in the course of the proof. By 



Proposition ! 1.81 there exists Cq > such that 



Since w s j[t=Q 
(7.3) 



SUp ||lli(t)|| E 3 + Sup ||ti 2 (t)|| E 3 < Co- 

te[o,T] te[o,T] 
and e C([0, r e ]; S 3 ) for some t £ , we can find t £ > such that 

K(*)IIe» + K(*)IIe» <Co 



for ^ t t £ . So long as d7.3l ) holds, we perform energy estimates, to show that d6.8b is 
true, with a constant C independent of e. It will follow that up to choosing e G (0, So] with 
eo > sufficiently small, d7.3l ) holds for t € [0, T], which yields Theorem |6.1| 



Notation. For two positive numbers a e and 6 e , the notation a £ < b £ means that there exists 
C > independent of e such that for all e € (0, 1], a e ^ C6 e . 

Note that so long as d7.31 l holds, similarly to the case of I s , Proposition ^. ll implies 



WHt) 



< F 3/2 



W 3 >< 



hence 



2 Re Wf ) ii 



<V~e. 



Therefore, the last line in (17.2b . viewed as a source term, is O(^fe) in E 3 , so long as ( 17.3b 
holds. The other terms in (17.2b can then be considered as linear terms, in view of the 
application of Gronwall's Lemma. 
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We write 

(*J,diag * l"j| 2 ) "J ~ (Kl diag * \ Uj \ 2 ) Uj = (Kl diag * \u £ \ 2 ) w £ 

+ feag * (K-l 2 - Kl 2 )) % 
= (^,dia g *l^| 2 )^ 

2ReujWj)) Uj, 



(«i,diag*(K| 2 



and a similar relation for the off-diagonal kernel. We develop the general convolution, 
where K £ is of the form (16.4b : 

[K £ *f)9 = (JJ (V-*, V 2 K (6(y - z)^e) (y - z)) (1 - G)d6f(z)dz) g 
= (y, (JJ^ (1 - 9)V 2 K (9(y - z)y/e) d6f(z)dz^j y\ g 
+ ^JJ (z, V 2 K (8(y - z)Je) z) (1 - 8)d9f{z)dz^ g 
- 2 (y, JJ^ (1 - 9)V 2 K (6(y - z)y/e) zdOf{z)dz^ g. 



The same computation is available for (16.5b . with heavier notations, so we leave it out. 
From this we readily infer 

IK^diag * (Kl 2 +2Ret^f)) % || s3 < \\w% 3 + ||^|| S 3 < H^He, 

where we have used Proposition 11.81 Cauchy-Schwarz inequality, and d7.31 > for the last 
estimate. We can infer an L 2 estimate for w £ : since all the terms of the form K £ * \u £ \ 2 
are real valued, the standard energy estimate yields 



K(*)IU a < / llfeag * (Kl 2 + 2Beu jW ?)) Uj \\ ds 
Jo 

+ / \\( K j,oS * (K+il 2 +2Rb%+i^+i)) u j\\ L 2 ds 
Jo 



1 -(2ReW £ ) 



L' 2 



ds+ / \\ P £ 3 (s)\\ L 2ds 



< 



(|K(S)|| S 3 + ||Wj + 1 (s)|| S 3) ds+ y/id 



To pass from this L 2 estimate to a S 3 estimate, we have to assess the action of the operators 
of multiplication by yj and V w on J7.2I ). First, V yj commutes with the left hand side of 
( 17.21 i. but not with the right hand side. We write 



• a w £ . 



I0I=3-|"I 



The first term vanishes in an L 2 estimate of dl em w £ , and in view of ( 16. 2K for all multi- 
indices a, /3 with s$ \a\ 2 and \/3\ = 3 - \a\, 

\\d' 3 V £ d a w £ \\ L2 < ||w||| E 3. 

Remark 7.1. The presence of the potential V £ , which is morally a time dependent harmonic 
potential, forces us to work in £ 3 , and not simply in H 3 (R d ): this is a standard feature 
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of such potentials, whose associated dynamics consists of rotations in phase space, so the 
regularity/decay of the functions must be the same in space and in frequency. 

The terms ^|. diag * \ u j\ 2 ^j w j ar, d (^J.off * l"j+il 2 ) w j are treated similarly, and 
produce a term of the form real xd%. im Wj, plus a term controlled in L 2 by ||to|||s3, so long 
as (E3 holds. 

On the other hand, the multiplication by yj commutes with the right hand side of ( 17.2l i. 
but not with the left hand side: 



id t + -A,y 



V, 



so the commutation errors for the equation satisfied by \yj \ 3 Wj consists of a linear combi- 
nation, with constant coefficients, of terms of the form y^d^Wj, with |a| + |/3| = 3. We 
end up with, so long as (17.3b holds: 



K(*)I| E 3< / {\\w%s)\\ s3 + \\w s j+1 (s)\\ i:3 )ds+ / V^ds. 
Jo Jo 

Gronwall's Lemma yields (16.8b . hence Theorem l6.ll 

8. Second order expansion and limiting phase shifts 

In view of Theorem 16.11 the phase shifts 8j are such that 9^ are bounded on [0, T], 
uniformly in e G (0, e a ], since \u e -\ 2 = \ii e -\ 2 = \uj\ 2 + 0(y/e). To study the limit of 0^ as 
e — > 0, we need to perform a second order expansion of (uf , u|) as e — > 0, to understand 
the contribution of order y/e. Therefore, we seek 

(8.1) u s j = Uj + *yiuf ] + 0(e). 

Remark 8.1. An error term of order O(e) is natural, since one could actually seek a more 
general asymptotic expansion to arbitrary order, of the form 

N 

(8.2) = Uj + ]T e^uf + O (e^' 2 ) . 

Resuming the arguments presented in Section [6] we see that formally, the last line in 
(16. 6t is 0(e°°), if we work with an infinite regularity. To get a second order approximation 
of u^, we simply need to prove that this term is 0(e), but we can certainly not perform the 
study with only an 0(^/e) information like we did in order to establish Theorem |6. 11 To 
compute the limit of 0|, we need to establish the asymptotic behavior of Uj up to 0{e) in 
S, and not only in L 2 , so we make an extra regularity assumption. We remark that if in 
Theorem l6.ll we require ai, E S 7 , with 

V e C* 7 (R+ x R d ; R), and df V € L°° (R + x R d ) , 2 < |/3| < 7, 
K e ^°°(R d ;R), 

then the conclusions of Theorem |6.1| can be replaced by: (uf , tt|) £ C([0, T]; E 4 ) 2 and 

sup ||uf (t) - ui(t)|| S 4 + sup \\u\{t) - ii 2 (t)|| S 4 < Cyfe. 
te[o,T] te[o,T] 

In particular, (wf,?/!) € C7([0, T]; E 4 ) 2 uniformly for e £ (0, £q]. Thanks to Proposi- 
tion l3.ll this enables us to claim that in ( 16.71 ). 



Wf(t) 
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hence 



2 Re Wf u 



To derive an equation for the corrector Uj, we plug d8.1| l into d6.7t , and discard all the 
terms which are, at least formally, 0(e), including thus the last line. The term correspond- 
ing to the power -Je yields: 



< 



•fl (!) 



-AtiS" = VW + (if° diag * K| 2 ) vf> + (K° oS * | Uj+1 | 2 ) 



(8.3) 



+ 2 (tf° diag * Re (ujuty) Uj + 2 (tf° off * Re Uj 

+ + (^J,diag * |Uj| 2 ) Uj + (/Cj.off * |u J + l| 2 ) Uj, 



with Cauchy data w^_ = 0, where we have denoted the third order Taylor expansions 

Vi(t,y) = ^ 3 V{t, qj {t))yyy, 

%diag(2/) = ^V 3 X(0)y • y • y, 
1 

^Cj,off(t, y) = K (qj(t) - q j+ i(t)) y-y-y. 

These equations are naturally linear in the unknown (ttj , *4 ). ^ n view of Proposi- 
tion [O] the last line in ( 18.31 . which corresponds to a source term, belongs to C(R+; E 4 ). 
This non-trivial source term makes Uj non-zero. Even though (18.3b is a linear system, it 
seems easier to prove that it has a unique solution, by proceeding in the same way as in the 
proof of Proposition ll.8l (see Section|5]l. We have: 

Proposition 8.2. Suppose that ax, <i2 £ E 7 , and 

V G C 7 (R+ x R d ; R), and d?V G L°° (R+ x R d ) , 2 sC \/3\ < 7, 
if G W 7 '°°(R d ;R). 

77ze« (18 - 3b has a unique solution (u^\ u^j £ C(R+; S 4 ). 

Denote by u| = Uj + ^/eUj the second order approximate solution, and by = 
Uj — Vj the corresponding error term. It satisfies Wj, t _ = 0, and 



id t w* + -AtSf = Vfi3] + (iq diag * |^| 2 ) u] - (iq diag * |^| 2 ) v] 
+ (K? off * |^ +1 | 2 ) u) - (Kl oS * |^ +1 | 2 ) «j 



where the new source term is such that 



sup ||p|(i)|| s < Ce. 
te[o,T] 



Resuming the energy estimates used in Section [7] we infer: 
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Proposition 8.3. Let T > 0. Under the assumptions of Proposition^^ there exists eq > 
such that for e G (0, Eq], ( 16.71 ) has a unique solution [u\ , uf) € C([0, T]; £ 4 ) 2 . Moreover, 
there exists C independent of e € (0, £q] smc/i f/iaf 



sup 

te[o,T] 



(f(t)-«i(t)- Veui X) (*) 



+ sup 

E te[o,T] 



«a(*)-«2(«)-Vi«4 1) (<) 



^ Ce. 



Note that unlike in the proof of Theorem 16. II no bootstrap argument is needed at this 
stage, since we already have uniform estimates for u e - , Uj , Uj in C( [0, T] ; E 4 ). We readily 
infer: 

^(t) = ^(t) + 0(Vi) inL^dCT]), 

where 0j is given by 



(8.4) 



VAT(O) • ( 2 Re / zu j {s,z)uf\s,z)dz ) ds 



+ y Vif(gj(s) — qj + x(s)) ■ ^2Re ^ zuj+i(s, z)u^l x {s, z)dz^j ds. 



We have obviously 9j G C x ([0,T]), and 0j(O) = 0^(0) = 0. To see that 9j G C 2 ([0,T]), 
in view of the Cauchy-Schwarz inequality, and since Uj, Uj G C([0, T]; S), it suffices 
that verify that Uj,u^ G C 1 ([0, T]; L 2 ). This property is a straightforward consequence 



of Equations (I1.13l l and ( 18.3b . in view of the regularity of Uj and uj 1 ' . This completes the 



proof of Theorem 1 1.9 1 

To conclude, we check that the phase shifts 9j are non-trivial in general, by computing 
their initial second order derivatives: since u }} „ = 0, 

]\t=0 

Sj(0) = V-ftT(O) • ^2 Re J za j (z)d t uf ) (0,z)dz 

+ VK{ qj {0)-q j+1 (0))- ^2Re j za j+ i{z)d t ufl x {Q,z)dz 

From (18. 3I >, we have 

id t uf\Q, y) = (V 3 (0,y) + /C,-,di ag * | « j | 2 +/C il0ff * |a j+ i| 2 ) a,,(y) 
so the first line in the expression for 9j (0) is zero, and 

§ j (0)=VK(q j (p)-q j+ i(p)) 

+ VK( qj (0)-q j+1 (0)) 

+ VK( qj (0)-q j+1 (0)) 
Therefore in general 9j ^ 0. 

Remark 8.4 (Instability). The fact that it is necessary to analyze an 0(y/e) correction to 
(u x , 1*2) to compute (9\ , #2) implies the existence of instabilities at the semi-classical level. 
Typically, a perturbation of the initial data at order e 7 with < 7 < 1/2 will affect the 
leading order behavior of u £ in L 2 (for the strong topology) for some time < t £ — > 0. 
On the other hand, since the 0j's are purely time dependent, the Wigner measure are not 



zVj(0, z)2 Im (aj +x aj) {z)dz 
(£j,dia g * Kl 2 ) 2Im(a J+ iOj) (z)dz^j 
(fcj,oS * K+i| 2 ) 2lm(a j+ iaj) (z)dz 
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affected by this phenomenon. Since the approach to describe this instability is the same as 
in iflOl . we simply refer to that paper for more details. 

Acknowledgement. I am grateful to Alexander Mielke for pointing out Remark [L7l 
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